arXiv: 1507.04341vl [math.PR] 15 Jul 2015 


Activated Random Walks* 


Leonardo T. Rolla 


Contents 

1 Overview [2] 

1.1 The Activated Random Walk reaction-diffusion model. |2] 

1.2 Infinite conservative system, fixation, phase transition. I3] 

1.3 Physical motivation: “self-organized criticality” . E] 

1.4 Predictions . 0] 

1.5 Difficulties. E] 

1.6 Results. E] 

2 Definitions and Diaconis-Fulton construction [8] 

2.1 Notation. E] 

2.2 Diaconis-Fulton construction. . M 

3 One-dimensional counting arguments 1131 

4 Exploring the instructions in advance 1161 

5 Higher-dimensional arguments 1191 

6 A multi-scale argument [24l 

7 Arguments using particle-wise constructions 1271 

7.1 Preliminaries . EZl 

7.2 Results. EHl 

References 1321 


*Minicourse given at the workshop “Activated Random Walks, DLA, and related topics” at IMeRA- 
Marseille, March 2015. Preliminary version. July 15, 2015. 


1 













1 


Overview 


1.1 The Activated Random Walk reaction-diffusion model 

The Activated Random Walk model is defined as follows. 

Particles sitting on the graph can be in state A for active or S for passive. Each 
active particle, that is, each particle in the A state, performs a continuous-time random 
walk with jump rate Da = 1 and with translation-invariant jump distribution given by a 
probability p(-) on 7/. 

Several active particles can be at the same site, and they do not interact among themselves. 
When a particle is alone, it may become passive, a transition denoted by A —)• S', which 
occurs at a sleeping rate 0 < A ^ cx). 

In other words, each particle carries two clocks, one for jumping and one for sleeping. Once 
a particle is passive, it stops moving, i.e., it has jump rate Dg = 0, and it remains passive 
until the instant when another particle is present at the same vertex. At such an instant 
the particle which is in S state flips to the A state, giving the transition A -|- S' —)■ 2A. 

If the clock rings for a particle to sleep while it shares a vertex with other particles, the 
tentative transformation A —)• S' is overridden by the instantaneous reaction A-|- S' —)■ 2A, 
so this attempt to sleep has no effect on the system state. 


Jump distribution: x ^ x + z with probability p{z). 

Diffusion: jump at rates Da = 1 and Dg = 0. No interaction. 

Reactions: A —)■ S' at rate 0<A^cx), Aj-S'—)• 2A at rate oo. 

Notation: rjt{.x) denotes both the number and type of particles at site x at time t. 


A particle in the S' state stands still forever if no other particle ever visits the vertex 
where it is located. At the extreme case A = J-cx), when a particle visits an empty site, 
it becomes passive instantaneously. This case is equivalent to internal diffusion-limited 
aggregation with infinitely many sources. 

We have described local rules for the system to evolve. In order to fully describe the 
system, we need to specify on which subset of Z*^ this dynamics will occur, what are the 
boundary conditions, and the initial state at t = 0. 
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1.2 Infinite conservative system, fixation, phase transition 

Consider a system running on the whole graph and such that the initial conhguration 
rjQ is i.i.d. Poisson with parameter fi. 

A phase transition in this systems arises from a conflict between a spread of activity and 
a tendency for this activity to die out, and the transition point separates an active and 
an absorbing phase in which the dynamics gets eventually extinct in any hnite region. 


We say that the system locally fixates if rit{x) is eventually constant for each x, oth¬ 
erwise we say that the system stays active. 


It can be shown using ergodicity that the probability that the system hxates is either 0 
or 1. Moreover, 


If the system hxates for given A and p,, then it hxates for larger A and smaller p. 


As a consequence, 


P(hxation) 


If, p <i Pc 

1^0, p ^ Pc 


where the critical density pc = Pc{^) G [0, cx)] is non-decreasing in A G [0, cx) . 


1.3 Physical motivation: “self-organized criticality” 

In contrast with the above description, consider now the following evolution. 

Driven-dissipative system. For large L, let the system run on the hnite box V = [0, 

New particles are added to the bulk of V at constant rate. When a particle hits the 
boundary of V it is killed. The reaction-dihusion dynamics is run at a much faster speed, 
so that the whole box is stabilized between two arrivals. We then let t —)■ cx) to reach 
stationarity, then take L 3> 1 to have a state on 

The relation between self-organized and ordinary criticality is understood as follows. 

On the one hand, “self-organized criticality” appears in the parameter-free, driven- 
dissipative evolution described above. In this dynamics, when the average density p 
inside the box is too small, mass tends to accumulate. When it is too large, there is 
intense activity and a substantial number of particles is absorbed at the boundary. With 
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this carefully designed mechanism, the model is attracted to a critieal state with an av¬ 
erage density given by 0 < /Xc < oo, though it was not explicitly tuned to this critical 
value. 

On the other hand, the corresponding conservative system in inhnite volume exhibits 
ordinary criticality in the sense that its dynamics fixate for /i < pc and do not fixate 
for n > idc, and moreover the critieal exponents of the finite-volume addition-relaxation 
dynamics are related to those of the conservative dynamics in infinite volume. 

1.4 Predictions 

The behavior of the ARW is expected to be the following. 

The critical density satisfies /ic —t 0 as A —)■ 0 and /Xc —t 1 as A —)■ cxo. See Figure [T] The 
value of fic should not depend on the particular /x-parametrized distribution of the initial 
configuration (geometric, Poisson, etc.). 

At /X = /Xc, the density of active particles vanishes as f —)■ cxd, but we conjecture that the 
system does not fixate in this case. 

The asymptotic decay of density of activity as f 3> 1 when /x = /Xc should obey a power 
law. Also, for the stationary regime, i.e., letting f —)■ cx) first, the density of activity 
should decay with a power law as 0 < /x — /Xc -C 1. More predictions can be found at 
jPRSin] . 



Figure 1: Predictions 


4 










1.5 Difficulties 


For the ARW, we would like to describe the critical behavior, the scaling relations and 
critical exponents, and whether the critical density is the same as the long-time limit 
attained in the driven-dissipative version. These questions are however far beyond the 
reach of current techniques. 

The hrst apparent difficulty lies in the fact that this system is not attractive (i.e., its 
evolution does not preserve monotonicity of conhgurations). This is overcome by con¬ 
sidering a Diaconis-Fulton kind of construction, or an explicit construction in terms of a 
collection of random walks, rather than the classical Harris graphical construction. These 
constructions allow for different kinds of arguments which have proven to be very useful. 

However, another particularity of this model brings severe difficulties, which is particle 
conservation. In particular, this seems to rule out any energy-entropy kind of argument. 
These arguments typically go as follows. One hrst identihes some structure that is intrinsic 
to the occurrence of events that conjecturally should not occur. The number of possible 
structures is then estimated, and when it is overwhelmed by the high energy cost needed to 
construct them, one can show that such events have vanishing probability. This approach 
has proven successful in perturbative statistical mechanics. However, the conservation of 
particles in the system gives rise to intricate long-range effects, which makes it difficult 
to hnd a suitable structure behind self-sustained activity for large periods of time. 


Lack of attractiv e n e ss . Overcome by using other constructions rather than Harris’. 
Conservation of particles. Rules out any “energy x entropy” kind of approach. 


1.6 Results 

We end this overview by recalling the existing results in the literature about the phase 
transition of the ARW on the inhnite lattice In the next sections we should go through 
the proof of all the results mentioned here. 

The results will be presented in the same framework, by making the following assumptions. 


We shall assume that: (i) the initial distribution is i.i.d. Poisson parametrized by p, 
and (ii) jumps are to nearest neighbors, that is, p{z) = 0 unless || 2 :|| = 1. 
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For each result, the assumption of i.i.d. initial condition and a Poisson law for its marginals 
can be more or less relevant depending on the proof. We cannot make any other assump¬ 
tions about the jump distribution, since some of the existing results have been proved 
only for biased jumps, whereas some others only in the complementary case. 

Each result can be generalized in different directions. The reader interested in minimal 
hypotheses is referred to the original articles. 

One-dimensional directed walks 

For the particular situation when particles only jump to their right, all the above predic¬ 
tions about the phase transition curve can be proved. More precisely, 

A 

and there is not hxation at /x = This appeared in |CRS14j although it was known 
before. Yet, much remains to be understood about its critical behavior. In particular, the 
scaling limit of the flow process has been studied in |CRS14j for the extreme case A = cx), 
but remains open elsewhere on the critical curve. 



Figure 2: Results for d = 1, totally asymmetric jumps 


One-dimensional systems 

In studying a different model, arguments found in |KS06j imply that for high enough /i 
the ARW with A = cxo model stays active. In terms of the phase diagram, fic < oo for 
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any A ^ oo. It was shown in |RS12] that there is no hxation at /r = 1, so in particnlar 
/ic ^ 1 for any A ^ cx). 

Finally, in |Tagl4| it was shown that when the jnmp distribntions is biased, pc < 1 for 
any A < cx), and pc —t 0 as A —?■ 0. It should be noted that the lower bound obtained 
for yUc tends to 1 if the bias is small, for any A > 0, and that the case of unbiased jumps 
remains open. 

The hrst result on hxation in this context appeared in |RS12] , showing that 

A 


he ^ 


1 + A' 


In terms of the phase diagram, he > 0 for all A, and /Xc —t 1 as A —)■ cx). 



Two and more dimensions 

In |KS06j . the fact that fic < oo was proved for any dimension. It was then shown that 
/ic ^ 1 by |Sheinj and |AGGinj . using different methods. This was extended in |GRS14] . 
who showed that there is no hxation at /x = 1. 

In |Tagl4| , it was shown that he < 1 for sufficiently biased jump distributions and some 
A > 0. Unlike the one-dimensional case, the proof gives neither pc < 1 for all A < oo, nor 
Pc —t 0 as A —)■ 0. The case of small or zero bias remains open. 

The hrst result on hxation for arbitrary dimensions appeared in [ShelOj . who shows that 
Pc > 0 when X = oo. This was extended in |CRS14] who show that Pc = 1 when A = cx). 
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Finally, in |ST14j . it was shown that when the jump distribution is unbiased, /Xc > 0 for 
any A > 0. The methods of |ST14j do not allow to prove that pc —t 1 as A —)• cx), though. 
The case of biased jumps and X < oo remains widely open. 



2 Definitions and Diaconis-Fulton construction 

2.1 Notation 

Let No = {0,1, 2,... } and = Nq U {p}. The state of the ARW at time t ^ 0 is given 
hy rjt E T, = In this setting, rit{x) = n means that, at time t, site x contains n 

active particles if n G No, or one passive particle if n = g. 

We turn N^ into an ordered set by letting 0<p<l<2<---. We also let |p| = 1, so 
\rit{x)\ counts the number of particles regardless of their state. We dehne |n] to be n if 
n ^ 1, and 0 if n = 0 or p, so |? 7 f(x)] counts the number of active particles. 

To add a particle to a site we dehne p+1 = 2, which represents the A + S' —)■ 2A transition. 
We also dehne 1 ■ g = g and n ■ g = n for n ^ 2, which represent the transitions A —)■ S' 
and 2A —)■ A + S' —2A, respectively. 

It will be convenient to consider as acceptable the operations p — 1 = 0 and g ■ g = g, 
although they do not appear in the dynamics. The operations 0 ■ g and 0 — 1, on the other 
hand, are not acceptable. 

The process evolves as follows. For each site x, a Poisson clock rings at rate (1 +A)|pt(a;)]. 
When this clock rings, the system goes through the transition p -E with probability 
















otherwise rj —)■ TxyTj with probability p{y — a;) The transitions are given by 


TxyTliz) = { 


r]{x) - 1 , 

v{y) + 1 , 

V{z), 


Z = X 

z = y 
otherwise, 


and Tx(>rj{z) 


{ g ■ r]{x), z = X 
r]{z), otherwise. 


We assume that Pq^x) G Nq for all x a.s., and use denote the law of {rit)t^o, where z/ 
denotes the distribution of rjo. 


2.2 Diaconis-Fulton construction 

The Diaconis-Fulton representation enables us to exploit the combinatorial nature of 
hxation. Due to particle exchangeability, this representation extracts precisely the part 
of the randomness that is relevant for the phase transition, focusing on the total number 
of jumps and leaving aside the order in which they take place. It is suitable for studying 
path traces, total occupation times, and hnal particle positions, but precludes the analysis 
of quantities for which the order of the jumps does matter, such as correlation functions 
or local shape properties. 

We summarize the main points of this section for later reference. 


- Toppling x is an operator that reduces the state at x and increases elsewhere. 

- Legal topplings are those actually performed by the system. 

- The order of topplings is irrelevant and we are free to choose it. 

~ Acceptable toppling is an artihcial operation that extends the dehnition of legal 
topplings, and provides upper bounds for activity. 

- On the other hand, legal topplings provide lower bounds for activity. 

~ This construction determines hxation or non-hxation for the continuous-time model. 

- There is a zero-one law for hxation. 


In this section, p denotes a conhguration, and we do not deal with a continuous-time 
evolution anymore. We say that site x is unstable for the configuration p if p{x) ^ 1. 
Unstable sites can topple, by applying Txy or Txq to p. Toppling an unstable site is legal. 
If p{x) = g, we still say that toppling x is acceptable. A legal toppling is also acceptable. 

Consider a held of instructions X = {T^’^)x£Z‘i,j£N- Later on we will choose X random, but 
for now X denotes a held that is hxed, and p denotes any conhguration. 
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Let h = {h{x)\x G count the number of topplings at each site, usually started from 
h = 0. The toppling operation at x is dehned by 

and is a short for 0)- Again, ^x is legal for rj if ri{x) ^ 1, and acceptable for rj 
if ri{x) ^ g. 

Sequences of topplings and local properties 

Let a = (xi,... ,Xk) denote a hnite sequence of sites. We dehne = ^Xk^x^_i ■ ■ ■ 

We say that is legal (respectively acceptable) for rj if $ 3 ,; is legal (respectively accept¬ 
able) for for each I = 1,..., k. In this case we say that a is a legal seguence 

(respectively acceptable seguence) of topplings for rj. 

Let TUa = (ma{x)-,x E be given by ma{x) = J2i ^xi=x-i the number of times the site x 
appears in a. We write ma ^ rng if ma{x) ^ rng{x) V x, and f) ^ ?7 if fj{x) ^ rj{x) V x. 
We also write (f), h) ^ ( 77 , h) ii fj ^ rj and h = h. 

Let a; be a site in and 77 , fj be conhgurations. 


1) Local Abelianness. If a and P are acceptable sequences of topplings for the 
conhguration rj, such that rUa = mg, then = ^/sV- 

2) Mass comes from outside. If a and P are acceptable sequences of topplings for 
rj such that ma{x) ^ mg{x) and ma{z) ^ mg{z) for z x, then ^ap{,x) ^ $^77(0;). 

3) Monotonicity of stability. If site x is unstable for the conhguration 77, and if 
fi{x) ^ ri{x), then x is also unstable for the conhguration fj. 

4 ) Monotonicity of topplings. If f) ^ 77 and <ha; is legal for 77, then <ha; is legal for 
fj and ^ ^xV- 


Proof. The two last properties are immediate. Suppose that ma = mg. 

For convenience, dehne the operators 77 © = 77 + 1 on N^,, we well as 770 = 77 — 1 and 
770 = 77 ■ ^7 on Np \ {0}. With this notation, whenever 77 © is acceptable (i.e., 77 7 ^ 0) 
we have 77 © © = 77 © 0. Analogously, whenever 77 © is acceptable (i.e., 77 7 ^ 0) we have 
77 © © = 77 © ©. Therefore, within any acceptable sequence of operations, replacing 0© 
and 0© by 0© and 0© yields an acceptable sequence with the same hnal outcome. 

Notice that ^aV{x) is given by ri{x) followed by a number of ©, © and ©’s. The number 
of times each operator appears is determined by X and m^ and is thus the same for 
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^I 3 ri{x). Their actual order depends on the sequence, but the internal order of the 0 and 
0 operators is determined by and is the same for both As a 

consequence, we can apply the above identities to move the 0 ’s to the left, yielding then 
identical sequences for $^ 77 ( 0 ;), which proves the hrst property. 

The second property follows from a similar argument. Suppose ma{x) ^ rriis^x) and 
^a{.z) ^ mi 3 {z) for z ^ X. Again is given by rj{x) followed by a number of 0, 0 

and 0 ’s. The number of times that operator 0 appears depends on ma{z), z ^ x, and is 
bigger than for <h^? 7 (a;), and the number of times that operators 0 and 0 appear depend 
on ma{x), and is smaller than for $^ 77 ( 0 ;). Pushing the 0 ’s to the left we get that 
is written in the same way as $^ 77 ( 0 ;), perhaps with a few extra 0 ’s in the beginning, and 
a few missing 0 and 0 ’s in the end, so $ 077 ( 0 ;) ^ $^ 77 ( 0 ;). □ 

Global properties 

What follows is valid for any model satisfying the above four properties. 

Let V denote a hnite subset of We say that 77 is stable in V if every a; G is stable 
for 77 . We say that a is contained in V, and write a C if every x appearing in a is an 
element of V. We say that a stabilizes 77 in if a is acceptable for 77 and $oh is stable 
in V. 


Least Action Principle. If a is an acceptable sequence of topplings that stabilizes 
77 in V, and C 1/ is a sequence of topplings that is legal for 77 , then ^ ttIq. 


Proof. Let /d C 1/ be legal and uia ^ tU/?. Write ft = {xi,..., Xk) and = (^xi ,..., Xj) 
for j ^ k. Let i = max{j : m^u) ^ ttIq,} < k and y = G V. Since /3 is legal, y is 
unstable in $^(o? 7 - But rupa) ^ tUq and mij(e){y) = maiy). By the Properties 2 and 3, y 
is unstable for $^77 and therefore a does not stabilize 77 in P. □ 

For each hnite set V and conhguration 77 , we dehne 

mvx = sup 771/3. 

ycV legal 

The Least Action Principle says that 

771/3 ^ ^V,rj ^ 

for any legal sequence fd contained in V and any acceptable sequence a stabilizing 77 in 
V. This provides a very good source of lower and upper bounds for mv^- 
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Global Abelianness. If a and /3 are both legal toppling seqnences for rj that are 
contained in V and stabilize r] in V, then ma = rriis = mv,rj- In particnlar, = ^pV- 


Proof. Apply the Least Action Principle in two directions: rrijs ^ ^ rua ^ frijs. □ 

We say that rj is stabilizable in V if there exists an acceptable seqnence a that stabilizes 
7] in V. Notice that this provides hnite npper bonnds for legal seqnences /? C 1/, which in 
tnrn implies the existence of a sequence contained in V that is both legal and stabilizing. 
Indeed, if one tries to perform legal topplings in V indehnitely, on the one hand one has 
to eventually stop since there is a hnite upper bound, and on the other hand this can only 
be stopped if there are no more unstable sites. 


Monotonicity, li V C V and rj ^ rj, then mv^n ^ 


Proof. Let /3 C 1/ be legal for rj. By successively applying Properties 3 and 4, /3 is also 
legal for fj. Since (3 <PV in this case, the inequality follows from the dehnition of mv,r)- □ 

By monotonicity, the limit 

rUj, = lim 

exists and does not depend on the particular sequence 1/ f (the limit is given by the 
supremum over hnite V). A conhguration rj is said to be stabilizable if mri{x) < oo for 
every x & 7/. 

Fixation for the stochastic dynamics 

Recall that denotes the law of the process {rit)t^o with values on Assume that 
distribution n of po on hlo"* ergodic and has hnite density u^rj^O)) < oo. We further 
assume that the support of p(-) generates and not a sublattice. 

To avoid extra notation, we dehne the held X on the same probability space P*^, indepen¬ 
dently of rjo, and distributed as follows. For each x G and j G N, choose 
with probability or r^g with probability independently over x and j. 

Stochastic Dynamics. P'^(hxation of {r]t)t^o) = < oo) = 0 or 1. 


For the proof, the reader is referred to |RS12l v2 on arXiv]. 
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3 One-dimensional counting arguments 


|CRS14] {d = 1. Directed walks). fXc = No fixation for jj, = /ic- 


Proof. Fixation is equivalent to m^Q(O) < oo almost surely. This is equivalent to mvj^,rioW 
being tight as L —)■ cxo, where Vl = {—L ,..., 0}. The latter is equivalent to tightness of 
the number of particles which jump into 0 when Vl is stabilized with legal topplings. See 
Figure 
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Figure 5: The set Vl, part of the conhguration t], and arrow indicating the unknown 
number of particles which jump into 0 when stabilizing [—L, 0] 


Let us stabilize Vl by exhausting each site, from left to right. More precisely, topple site 
X = —L until each of the rio{—L) particles either moves to x = —L + 1 or sleeps, and let Yq 
be indicator of the event that the last particle remained sleeping on x = —L. Conditioned 
on TjQ^—L), the distribution of Yq is Bernoulli with parameter (in case rjo^—L) = 0, 
sample Iq independently of anything else). The number of particles which jump from 
X = —L to X = —L + 1 is given by Ni := [?7o(— L) — To]"'"- See Figure]^ 



Figure 6: Stabilizing [—L, 0] from left to right. Disks represent particles initially present 
at each site, and circles represent particles that arrive from the left. Red boxes indicate 
particles which became passive when they were left alone. 
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Note that, after stabilizing x = —L, there are iVi + r]o{—L + 1) particles at x = —L + 1. 
Let the site x = —L + 1 topple until it is stable, and denote by Yi the indicator of the 
event that the last particle remained passive on x = —L + 1. The number of particles 
which jump from x = —L + 1 to x = —L + 2 is given by N 2 '■= [iVi + r]Q{—L + 1) — Ti]^. 
By iterating this procedure, the number Nj+i of particles which jump from x = —L + i 
into X = —L + i + 1 after stabilizing x = —L, —L + 1,..., —L + z is given by 

where Nq = 0. The number of particles which jump into 0 while stabilizing Vl equals Nl. 

Note that the sequence (iVi)i=o,i,...,L is a random walk on No, with independent increments 
distributed as rjo^x) — Y, reflected at 0. So the relevant quantity is 

E[po(-T + fc) - Yfc] = /i - 

If h < fhs walk is positive recurrent, which implies tightness of iV*, and thus hxation. 
If p > the walk is transient, and P[iVi ^ f (/i — is large for large L, so there 

is no hxation. If /i = the walk is null-recurrent, so there is no tightness and thus no 
hxation. □ 


{d = 1). No hxation when p = 1. 


Proof. Let p = 1. By the CLT, the probability that rjo contains at least L + 2'/L particles 
on Vl = [0, L] is at least 26 > 0, uniformly on L. On this event, at least 2^/L particles 
will visit either x = 0 or x = L when we stabilize [0, L]. Therefore, 

26 ^ F(mv^,r^o{0) ^ +F(mv^,r,o{L) ^ V^ 2P(m^g(0) ^ Vl). 

Since this is true for any L, there is no hxation. □ 


|Tagl4| (d = 1, Biased walks). For any A < cx), pc < 1. Moreover, limA-s-o PcW = 0. 


Proof. Let F{X) be the probability that a walk which jumps at rate 1 and sleeps at rate 
A, but only sleeps at (—cx),0], is able to reach -|-cx) without sleeping. A more explicit 
dehnition is F{X) = E[(j^)^], where T counts how many steps a discrete-time random 
walk starting at x = 0 is found on (—cxo, 0]. We will show that pc ^ 1 — F{X). 
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Suppose /i > 1 — F{X). We will perform a sequence of legal topplings on [—2L, 0] as 
follows. 

In the first stage, we move each particle initially found on [—L, —1] until it is alone, in 
order to have a configuration on {0,1}^. If particles reach {0} they are stopped, but on 
[—2L, —L] we still move them. With high probability, none of the biased random walks 
will go further to the left than log^ L. Therefore, after this step at most log^ L particles 


—L — \og^L —L + 1—L + 2 —1 6 1 


















-L- 

log"L 


- 

L -L 

+ 1 -L 

+ 2 




- 

^ i 

i - 

L 


Figure 7: The two stages in the toppling procedure. After the second stage, a number of 
particles is found at 0. 
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will be outside [—L, 0 ]. 

We now move to the second stage. Let Nq denote the number of particles on Aq = [—L, 0] 
after this first step. Notice that 

At this point, all sites in this set, except site 0, have either 0 or 1 particle (*) 

Let X = —L and Ai = [—L + 1,0]. If x has no particles, then at Ai there are Ni = No 
particles, and Ai satishes (*). If x has one particle, move this particle until it either sleeps 
on [—oo, x], or finds an empty site on [x + 1, —1], or else it reaches {0}. At the end of this 
step, the number of particles present at Ai will be either A"o or A"o — 1. The latter happens 
with probability at most 1 — F{X). In any case the conhguration on Ai satisfies (*). 

Move to X = —L + 1, define A 2 = [—L + 2, 0], again move the particle found at x (if any) 
until it reaches an empty site. Again the number N 2 of particles in A 2 after this step is 
W — 1 with probability at most 1 — F{X), otherwise is equal to Ni, and A2 satishes (*). 
Then move to x = —L + 2, dehne A 3 , N 3 , and so on. 

Finally, when x = —1 we have the set Al = {0} containing Nl particles. By applying 
the LLN, both to the number of particles initially found on [—L, 0], and to the number of 
particles lost on the second stage, we have that Nl > (/r — e)L — log^ L—[l — F{X)+e]L S> 1 
with high probability, hnishing the proof. □ 

Notice that the upper bound for jjLc degenerates as the jumps become less biased. 


4 Exploring the instructions in advance 


|RS12] {d = 1). For any A > 0, /ic ^ 


Sketch of the proof. The proof uses an algorithm that tries to stabilize all the particles 
initially present in po, following the instructions in X, with the aid of acceptable topplings. 

After describing the algorithm, we will show that, whenever it is successful, it implies that 
771^1,(0) = 0. We finally show that the algorithm is successful with positive probability if 
fi < By the 0-1 law, this implies almost sure fixation of the ARW model. 

1. Description of the algorithm. 

The algorithm consists in applying a trapping procedure to each particle. This procedure 
explores X until it identifies a suitable trap (given by a carefully chosen sleep instruction) 
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for the particle. To do that, it follows the path that the particle would perform if we 
always toppled the site it occupies, and stops when the trap has been chosen. In the 
absence of a suitable trap, the algorithm fails. 

Remark that some of the explored instructions are actually not going to be used by the 
corresponding particle by the time it settles at the trap, leaving some corrupted sites that 
could interfere with the statistics of the subsequent steps. 

If there is a particle at 0, we declare the procedure unsuccessful and stop. Otherwise, 
label the initial positions of the particles on Z by ■ • • ^ x _3 ^ x _2 ^ x_i < 0 < Xi ^ 
X 2 ^ ^ . Let Oo = 0. 

We now describe the trapping procedure for particle Xk, k ^ 1. Suppose that the hrst k — 1 
traps have been successfully set up at positions 0 < Oi < 02 < • • ■ < 0^-2 < cik-i < Xk-i- 
The set [0, 0 ^- 1 ] contains all traps and corrupted sites found so far. 

The settling procedure starts with an exploration. Starting at examine and follow the 
instructions in X one by one (whenever a sleep instruction is found, the next instruction 
at the same site is to be examined). Follow this exploration until reaching Ok-i- 

Next we set up the trap. During the fc-th exploration, we are sure to visit every x G = 
[ofc-i + l,Xfe — 1]. Moreover, the last instruction explored at each x G Rfc is a jump to 
the left, see Figures and For each x G Bk, the second last instruction may be a sleep 
instruction. If this is not the case for any x G Rfc, we declare the procedure unsuccessful 
and stop. Otherwise, let be the leftmost site at which the second last instruction 
explored was a sleep instruction, and call this second last instruction the k-th trap. 

Since the trap is a sleep instruction found immediately before the last instruction, which 
is a jump to the left, we are sure that the exploration path has not gone to the right of 
Ok after hnding the trap, so all the corrupted sites will be in [ok-i + 1,0^], see Figure 
This procedure can be carried on indehnitely, as long as all the steps are successful. 

II. If the algorithm is successful, then m^g(O) = 0. 

We will show that, following the instructions of X, rjo is stabilized in I 4 = [x-n^Xr] with 
hnitely many acceptable topplings, without toppling 0 . 

Let us hrst stabilize the particle that starts at xi. To that end, we successively topple 
the sites found by the hrst explorer, until it reaches the trap at oi. At this moment the 
particle will become passive, and the site ai will be stable. Notice that this is acceptable 
since the particle is following the same path that the explorer did, even if it will sometimes 
be passive. 

Notice that, after the last visit to Oi, the explorer did not go further to the right, so when 
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Figure 8: First exploration path. It starts at position xi of the first particle and stops 
when it reaches the origin. The horizontal axis represents the lattice, and above each site 
X there is a sequence of instructions The bold arrows indicate the last jump found 

at each site x G — 1], and the bold cross indicates a sleep instruction found just 

before the last jump, this being the leftmost such cross, whose location dehnes oi. 

settling the first particle we use all the instructions examined so far, except some lying 
in [oo + 1, Oi]. Therefore, the same procedure can be applied to the second particle, as it 
will hnd the same instructions that determined the second exploration path. 

Notice also that the hrst particle does not visit 0, and the second particle neither visits 
0 nor Oi, thus it is settled without activating the hrst particle. Following the same 
procedure, the fc-th particle is settled at a*,, without ever visiting {0, oi, 02 ,..., Ofe-i}, 
for all A; = 1,..., n. After settling the n hrst particles in Z+, we perform the analogous 
procedure for the hrst n particles in Z_. 

This means that rjQ can be stabilized in Vn with hnitely many acceptable topplings, not 
necessarily in Vn, and never toppling the origin. By the Least Action Principle, mv„^n{0) = 
0. Since it holds for all n G N and 14 t ^ ^is n —)■ cx), this gives m^(0) = 0. 



Figure 9: Second exploration path. It starts at position X 2 of the second particle and 
stops when it reaches Oi. The regions in gray indicate the instructions already examined 
by the hrst explorer. The dark gray contains instructions examined but not used, whose 
locations determine the set of corrupted sites. 
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III. The algorithm is successful with positive probability. 

For each site a; G i?i, the probability of Fueling a sleep instruction just before its last 
jump equals and this happens independently of the path and independently for each 
site. Thus, ai — oq is a random variable distributed as a geometric with parameter 
truncated at xi — qq. 

Since no corrupted sites were left outside [oq + 1, ai] in the hrst exploration, the inter¬ 
distance 02 — oi is independent of oi. Moreover, its distribution is also geometric with 
parameter The same holds for 03 — 02 , 04 — 03 , etc. By the law of large numbers, 
ttn ~ On the other hand, Xn ~ iJ,~^n. Therefore, if /x < the probability that 

Ofc < Xk for all k is positive. This event, in turn, implies success of the algorithm. □ 


5 Higher-dimensional arguments 

Let = (X^)n=o,i,2,... denote a random walk starting at z G and jumping according 
to p(-), the superindex is omitted when z = 0 . 

|ShelO] . fie ^ 1- 


Proof. Let Vn be the discrete ball of radius n around 0. Label and order all the particles 
initially found in Vn- We move the first particle until it is alone in a site, or until it exits 
the box. This can be done by subsequently toppling the site that contains the particle 
until one of these two conditions are met. If the particle is initially alone, there is nothing 
to do. After we are done with the hrst particle, we move the second particle until it is 
alone or exits the box. Then we move the third particle, and so on. 

The path performed by each particle in this procedure is distributed as a random walk 
truncated at the time when it reached an empty site (we do not care about what the 
particle does once it exits Vn). In order to complete the path, we sample a random walk 
starting from the point where the particle was stopped. This second part of the path, 
which is not actually performed by the particle, is called a ghost. See Figure [T 0 | 

Let W count the total number of walks that visit 0 before leaving W- Let L count 
the number of walks that visit 0 before leaving Vn, but which do so as ghosts. We are 
interested in counting the number of such walks that visit 0 as particles, which is given 
hjW-L. 

It may be complicated to determine the set of locations where ghosts start. But we can 
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use the fact that at most one ghost starts at each site of Vn- Let us start an artificial ghost 
from each site where no actual ghost has been started, as also illustrated in Figure [TOl so 
that exactly one ghost (artihcial or not) starts from each site. Let L ^ L denote total 
number of ghosts which visit 0 before leaving Vn- 



Figure 10: Particles which visit 0 before or after exiting Vn, ghosts started each time a 
particle hnd itself alone at a site, and artihcial ghosts started from other sites. 

Although W and L are far from being independent, we claim that each of them has 
variance bounded by a multiple of its expectation, and thus is concentrated around its 
mean. Finally, we also claim that 

E[hF] = /iE[Z] —)■ cx) as n —)■ cx). 

Hence, for /i > 1, the weak law of large numbers gives hF — L —>• oo in probability as 
n —)■ cxo, proving non-hxation. 

To complete the proof, let us prove the above elementary claims. First, 

m = x: where qx = visits 0 before exiting Vn\- 

x&V„ 

Also, E[hF|r7o] = J2xVo{^)Qx, so E[hF] = EE[LF|77o] = fM[L]. Write m = [n/S]. Estimat- 


20 







ing the above sum gives 

E[L] S ^ P[X'^ visits 0 before exiting + x] 

= ^ P[X° visits ^ before exiting Vm] 

zeVm 

= E [number of sites visited by before exiting Vm] —t oo 

as n —)■ oo. For the variance, hrst notice that L is given by the sum of independent indica¬ 
tor functions, so VL ^ EL. Conditioned on r^o, W is also given by the sum of independent 
indicator functions, so EV[hF|r7o] ^ EE[hF|77o] = EhF and VE[hF|77o] = J2x^[Vo{x)]qx = 
/i/i“^E[hF], where /r denotes hrst the variance and then the expectation of a Poisson. 
Finally, the law of total variance gives V[hF] ^ 2E[hF], hnishing the proof. □ 


rfagl4| (Biased walks). If the the walk is sufficiently directed, then there exist A > 0 
and /i < 1 for which there is no hxation. 


Proof. Denote the drift by n = J2zP{z)z e \ {0}. We assume that v ■ ei > 0 by 
symmetry, and that d = 2 for simplicity. We say that the walk is good if {Xf — z)-ei > 
0 for all n ^ 1. Let 

/C = P[(X„)„ is good] > 0. 

Notice that /C tends to 0 as the walk becomes less biased, and tends to 1 when the 
probability of jumping to the left is small. We will show that the system stays active if 



^(1 
i+aD 




/C > e-'^. 


This condition is meaningful when /C The result is weaker than in d = 1 because 

the argument requires some control on the position where particles exit a large box. 

Let Vn = {z = {x,y) G : a: G [—n, 1],?/ G [—Order these sites by writing 
Vn = {zi, Z 2 , ..., zn}, where N = 2n'^, observing that sites with smaller x must appear 
first. We now stabilize W as follows. 

Pick one particle initially present in zi, and move that particle until it either sleeps or exits 
Vn. If it sleeps in a given site of W, we leave that particle in that site and do not touch it 
anymore. In this case we start a ghost from that site so that the whole path associated to 
that particle is distributed as a random walk, part of which is performed by the particle, 

^For Bernoulli initial distribution, the result is a bit stronger. The term e~^ is replaced by 1 — /i. In 
this case, the above condition gives < 1 for any A < oo and any /C > 0 . 
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and part of it by the ghost. We say that the particle is good if the corresponding walk is 
good. 

Repeat the same procedure for all particles initially present in zi, except for the last 
particle. When a single particle is left at zi, its hrst action may be to sleep, otherwise it 
is to jump. If it jumps we treat it as the previous particles, so that we assign a walk to 
it (part of which may consist of a ghost walk). If it sleeps, we declare that particle dead, 
and do not even assign it a walk. 

After we are hnished with all particles initially present in zi, we repeat the same action 
with particles found on Z 2 , then z^, and so on. 

Sites which contain a dead particle are those initially containing at least one particle and 
such that the last particle sleeps instead of performing its hrst jump, so the density of 
dead particles is (1 — The total number of walks is given by the number of 

particles minus the number of dead particles, so it has density 

h :=h- 

When the particle is good, its walk is always to the right of its initial position, and the 
particles initially present on this region have not yet been moved. Therefore, a good 
particle can only sleep at the set of sites that were initially empty, which has density e~^. 

These observations explain the requirement fiK, > e~^. Let us start a ghost walk at every 
site that was initially empty, even if no particle was actually stopped there. This way we 
obtain an upper bound for the set of ghosts associated to good particles. It thus follows 
by the LLN that at least (/i/C — e~^)N — o{N) good particles exit 14 before sleeping. 
However, we want to estimate the number of particles that exit 14 through site 0. 

Let W count the number of good walks that exit 14 through 0 in this procedure. Let L 
count the number of good walks that exit 14 through 0, but which do so as ghosts. We are 
interested in the number W — L of good particles that exit 14 through 0 before sleeping. 
In principle it is hard to estimate L, so we consider the number L ^ L of ghosts that exit 
14 through 0 , when ghosts are started from each z ^ Vn that was initially empty. 

Since both W and L are given by sums of independent variables whose variance is com¬ 
parable to the expectation (see the proof of [ShelOj above), it suffices to show that 

E[I4] ^ flJCn — o{n) and E[L] 4 e~^n 

in order to have IT —L—)-cxoas?7,—)■ oo, in probability, which in turn implies non-hxation. 


22 




We use a simple argument of symmetry. First, writing Gk = {{x,y) ■. x = k}, 


— 1 

^[L] = X! X! exits 14 through 0] 

x=-ny=-rfi 
-1 

4 e“^ ^ ^ P[X*’^ reaches Go at 0] 

x=—ny£Z 
-1 

= e~'^ ^ ^ P[X'^’° reaches Gq at (0, —y)] = e~^n. 

x=—n ygZ 

Finally, writing Hk = {{x,y) ■. y = ±fc}, 

— 1 

E[1F] = E E pP[W^’^ is good and exits 14 through 0] 

x=-n y=-n^ 

— 1 

E E P[X'^’^ is good and reaches Gq at 0, before reaching y + 

x=-n y=-n'^ 

— 1 'n? 

= /^ E E P[X^’° is good and reaches Go at (0, —y), before reaching if„ 2 ] 

x=-n y=-rfi 
-1 

= /i ^ P[X^’° is good and reaches Go before Fr„ 2 ] 

x=—n 

-1 

4 /h E is good and reaches Gx+n before Fr„ 2 ] 

x=—n 

= llnF[X is good and reaches G„ before Ffn^]- 

Since ^ n as A; —>■ cx), the last probability tends to /C as n —>■ cx). This proves the 
desired inequalities, thus hnishing the proof. □ 


|Shein] . For X = oo, yc > 0- 


Proof. We will in fact show that, if y is small, then for almost every initial conhguration 
?7o, for any realization of X, only hnitely many particles can visit the origin. The idea is 
to assume that an adversary is trying to bring as many particles to the origin as possible, 
being able to move the particles to any nearest neighbor rather than following the random 
jumps. The proof uses the percolative structure of low-density sets. 

Given rjQ, for each hnite set V C dehne its weight by wiV) = Y.x&v ^ finite 

connected set 1/ C is internally fillable if w{V) 4 |ld|. 
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Consider the decomposition of 1/ into clusters C given by 

C(x) = y C : C is finite and connected, x G V, and w(V) ^ \V\ or \V\ = l|. 
These are called clusters because C{x) = C{y) whenever y G C{x). 

The main observation is that, regardless of X, for a particle starting at y to visit site x, it 
is necessary that some hnite connected set V containing both x and y is internally hllable. 
In other words, it is necessary that y G C(x). 

Therefore, in order to show hxation, it suffices to prove that, with positive probability, 
|C(x)| < cxo. This condition is equivalent to 

sup||I/| : ^ l| < cx). 

It follows from results in |Marn2] about “greedy lattice animals” that, for some C > 0, 

a.s. lim max ^ C'E[? 7 o( 0 )'^'*~^]. 

n^oov.\V\=n \V\ y J 

Since we are assuming Poisson initial distribution, the bound can be made less than 1 by 
choosing /i small enough, which finishes the proof. □ 


6 A multi-scale argument 


|ST14] (Unbiased walks). For any A > 0, /Xc > 0. 


Overview of the proof. The proof is too long and technical for these notes. We give an 
overview of the general strategy, omitting many delicate points. 

The main step is to show that an initial configuration restricted to a very large box 
stabilizes within a slightly larger box, which high probability. This is proved by recursion 
on the scale of the box, so the proof looks very little into the details of the actual ARW 
dynamics. In a sense, this kind of approach fits to our intuition that no matter how big 
a defect is, it will only affect a neighborhood of comparable size. 

The box at scale k is a, cube of length Lk, defined as follows. Let 7 = 1/10, Xq = and 


Notice that Lk increases as a doubly exponential of k. We also let Rk+i = [X/jXfc as 
an intermediate scale between Lk and Lk+i- In Figure 11 we see an inner box an 

intermediate box and a full box Bk+i of level k + 1. 
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^k + 1 


Figure 11: Boxes and scales; <C Rk+i ^ L^+i. 


Let pk denote the probability that, starting from a Poisson in some particle exits B^- 
That Pfc —)■ 0 fast as A: —?• oo follows from the recursion relation 

T 2d 

Pfc+l ^ Jjld Pk + ^k+l, 

consisting of a combinatorial term, the probability that stabilization fails twice at 
scale k, and the probability e^+i that something goes wrong at scale A: + 1. Indeed, if pk^ 
is small enough and Cfc —)■ 0 fast enough, then the square power above beats the 1 + 27 
power in the dehnition of Lk, and pk vanishes doubly-exponentially fast in k. 

Let us describe some aspects of this recursion step, depicted in Figure [T^ 

States in light gray have Poisson product distribution with the right density. They are 
restricted to the inner box of level A: + 1 for the initial conhguration, and the inner boxes 
of level k for the “sieved conhgurations”. States in dark gray are stable states, typically 
attained by the dynamics. States in gray are “balanced conhgurations”. Thick arrows 


Initial configuration Balanced configuration Sieved configuration 
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Failure at level k 


Balanced configuration 



Figure 12: Illustrative diagram of events for the recursion relation 


25 




























































































































represent typical events, while thin arrows represent events of low probability, either Ck+i 
or pk. 

Starting fresh. To let the dynamics rnn on boxes of the previous scale and use recursion, 
it is important to start with a Poisson product distribution within their inner boxes. This 
is achieved by a sieving procedure described below. 

Worst case scenario. If the dynamics fail to stabilize all the boxes of level k, the 
conhguration inside these boxes is no longer Poisson. In the absence of any useful knowl¬ 
edge about the resulting distribution of particles in this case, we use only the fact that 
the total number of particles within each box is still a Poisson random variable, and thus 
cannot be much larger than its mean. A balanced configuration is such that the number 
of particles within each box of level k is appropriately bounded. 

Sieving procedure. Starting from a balanced conhgurations, we move each particle for 
quite enough time until their relative position within the box of level k that contains it 
has mixed. If the particle happens not to be in its inner box, we repeat the procedure 
until it is. This reshuffling with sieving results in a state that can be coupled with an 
i.i.d. Poisson conhguration with high probability. This is one of the heaviest statements 
in |ST14] and not simple to prove. In order for this coupling to be possible, a slight 
increase in the density is necessary, analogous to the sprinkling technique in percolation 
(this increase should decay just fast enough so that it is summable over k). 

The chain of events. By hypothesis we start with a Poisson product measure inside the 
inner box The hrst good event is that such a conhguration is balanced, that is, each 

box Bk has a bounded number of particles. We then let these particles remix by forcing 
them to jump. During this procedure they can reach the intermediate box, bot not the full 
box. The second good event is that the resulting conhguration is properly sieved. We now 
let the evolution run normally within each box Bk, and the third good event is that each 
box stabilizes nicely without letting particles leave. It may happen that some of these 
boxes of level k is not stabilized. In this case, the conhguration is still balanced with high 
probability. We remix them again, now obtaining a sieved conhguration in the full box 
Bk+\. The system is given a second chance to stabilize, which should happen typically, 
but may fail again if again some of the level-fc boxes does not stabilize as expected. 

In the proof there are several aspects to keep under control, and several delicate statements 
that we omit here. The above description is not intended to let the reader understand 
the sketch of proof, but hopefully gives a general havor of the main argument. For all the 
details the reader is referred to the original article. □ 
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7 Arguments using particle-wise constructions 


The techniques presented so far used the Diaconis-Fulton construction and its properties, 
as described in |RS12] . In the particle-wise construction, the randomness of the jumps 
is not attached to the sites, but to the particles. In this section we consider techniques 
which use the particle-wise construction, or sometimes a combination of both. 

7.1 Preliminaries 

System with labeled particles Each existing particle at time f = 0 is assigned a label 
{x,j), where x & 1/ denotes its starting position and j = 0 ,... ,riQ{x) — 1 distinguish 
particles starting at the same site x. Let = {Y^ be given by the position of 
particle (x, j) at each time t. We write Y^’^ = T ^1/ \i rio{x) ^ j. Let 7 *’-^ = {t))t^o 

be given by 7 *’-^(t) = 1 if particle {x,j) is active at time t or p if it is passive. Write 
Y = and 7 = Then the triple = ( 70 , Y, 7) describes the whole 

evolution of the system. 

Whereas the process given by 

Vtiz) = 

^ j<Vo{x) 

only counts the number of particles at a given site at a given time, having each particle 
labeled gives a lot more information and allows different techniques to be employed. 

In a system whose initial conhguration rjo contains hnitely many particles, the evolution 
described above is always well dehned, since it is simply a continuous-time Markov chain 
on a countable state-space. Many different constructions will produce with the correct 
distribution. 

Particle-wise randomness Assign to each particle {x,j) a continuous-time walk 
independently of anything else, as well as a Poisson clock C M_,_ 
according to which the particle will try to sleep. X^'^ is the path of the particle param¬ 
eterized by its inner time, which may be slowed down with respect to the system time, 
depending on the interaction with other particles. (In the ARW model, when the particle 
is passive, the inner time halts until it is re-activated by another particle.) For this reason, 
X^'^ will be called the putative trajectory of particle {x,j). These random elements will be 
denoted by ^ = ( 7 o,X,P). Let a^'^{t) denote the inner time of particle {x,i) at instant 
t. Then {Yt'^)t^o is given by Yt'^ = X^Y^^y 
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Construction of the infinite system When rjo contains finitely many particles only, 
^ is determined by ^ in the obvious way, and it works for a.e. X and 7^. The simultaneous 
construction on the whole space, is done on the sequences of balls B{y,n), n & N, centered 
at each site y G This family of sequences is countable and translation-invariant. For 
each ball, consider the evolution of the system with initial conhguration r/o ■ ^B{y,n)- We 
say that this construction is well defined if (i) for each x and j, and for each T > 0, both 
and are the same in the systems {rjo ■ t) for all but hnitely many n, 

and (ii) the limiting process does not depend on y. 

The above dehnition is convenient for two reasons. First, the family of hnite approxima¬ 
tions to is countable, which makes it possible to prove that the construction is a.s. well 
dehned. Second, this family is translation invariant, so that, whenever the construction is 
well dehned, it is translation covariant. More precisely, for any translation 6 of Z^, if <^(^) 
is well dehned, then C(^^) is well dehned and equals 9C{^), where 9^ = {9rio,9'K,91^), 
etc. In particular, this implies that the system is ergodic, satishes the mass transport 
principle, and can be approximated by hnite systems regardless of which construction is 
used. 


The Mass-Transport Principle Let m : Z'^ x Z'^ —)• IR+ be a translation-invariant 
random function, that is, m{x,y,^) = m{9x,9y,9^) for any translation 9. The Mass 
Transport Principle is given by 


y 


E[^m{y,x) 

y 


Informally, the MTP says that the amount of mass transmitted from a vertex x is equal to 
the amount of mass entering x. It seems like an innocent and perhaps obvious identity, but 
its strength lies in its versatility, since it holds for any such function. The proof consists 
simply on re-indexing the sum and using translation invariance. See |LP| Chapter 8] for 
applications and generalizations to other settings. 


The particle-hole model Let us introduce a related model that will be useful in the 
sequel. Particles perform continuous-time random walks independently of each other. 
Sites not containing any particle are called holes. When a particle is alone at some site, 
it settles there forever, hlling the corresponding hole. After the hole has been hlled, the 
site becomes available for other particles to go through. If a site is occupied by several 
particles at t = 0“, we choose one of them uniformly to hll the hole at t = 0, and the 
other particles remain free to move. 






As in the ARW with A = oo, once a site has at least one particle, it will always retain one 
particle. The differences are (i) sites with n particles are toppled at rate n instead of n — 1 
and (ii) in the system with labeled particles, it is the hrst particle to arrive at a site that 
is retained, whereas in the ARW the particles can replace each other. Nevertheless, both 
models have the same Diaconis-Fulton representation, and therefore they are equivalent 
in terms of hxation at a; = 0. The Poisson clocks 7^ are not used in either system. 

7.2 Results 

Fixation as defined so far concerns the state of sites, and will be called site fixation. When 
each labeled particle eventually hxates we call that particle fixation. 


|AGG10j Assume that the particle-wise construction is a.s. well dehned. If some 
particles are not hxating, then sites are not hxating. Therefore, /ic ^ 1- 


Let us see how /ic ^ 1 follows from the main claim. First, it follows from the MTP that: 

The density of particles that fixate cannot he larger than 1. 

Indeed, let A{x,j,y) denote the event that particle {x,j) hxates at site y. Let m{x,y) = 
Jfj^A(x,j,y)- The density of particles that hxate is given by KJfym{x,y) and the density 
of sites eventually occupied by a hxating particle is given by KJfym{y,x). These are 
equal by the MTP. Since at most one particle can hxate at a given site, Jfyrn{y,x) ^ 1, 
proving the estimate. Now, if p > 1, there is a positive density of particles which do not 
hxate, which implies that there is no site hxation. Therefore, pc ^ 1. 

Proof. Let denote the event > j and {v,j) is non-hxating” and write A^ = A^'^. 

Suppose the probability that some particle is non-hxating is positive. Then there is j 
such that P(A°’-^) > 0. By interchangeability of particles, a := P(A°) = F{r]Q{0) > 
0, (0, 0) non-hxating) ^ F{t]o{0) > j, (0, 0) non-hxating) = P(A°’-^) > 0. 

As a warm up, notice that by the MTP the number Mt of non-hxating particles present at 
0 at time t satishes EiVA ^ cl, for any t > 0. Hence, the liminftEMj ^ a. If we had some 
control on YMt, we could conclude that hminf 4 P(Mt ^ 1) > 0, implying non-hxation. 
To control variance, we will use a local condition instead, plus extra randomization. 

Since the system is a measurable function of the randomness for any e > 0 there is 
fc G N such that the event A° can be e-approximated by some event that depends only 
on {riQ{x), ,V^) for ||x|| ^ k. Let B'^ denote the corresponding translation of the event 
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. When occurs, we say that particle (n, 0) is a candidate. It is a good candidate if 
also occurs, otherwise it is a bad candidate. 

Fix T > 0. For n G N, let us add more randomness to the system by sampling Z'" 
uniformly amongst the first n different sites in the putative trajectory X'”’^ after time T, 
independently over v. Define C{v,u) as the event and = m”. Let 


q{v,u) 


P(C'(r;,M) 


«) 


By the mass-transport principle. 


and Q{u) = q{v, u). 

V 


E[Q(n)] = ^P(C'(M,n)) = ^P(C(n,M)) = P(S^) =: 6 > 0. 

U U 


Notice that q{v,u) ^ K Notice also that q{v,u) and q{w,u) are independent if \w — u\ > 
2k. Using these two facts, it can be shown that, as n —)■ cxd, V[Q(r’)] —)■ 0, and thus 
Q{v) ^ b in probability. 

Let N{v) = J2u '^c(u,v) count the number of candidates for which = v. Then 


F{N{v) = 0\^)=l[{l-q{u,v)) 


^ e 


-Q(D 


—)■ e ^ in probability as n —)■ cx). 


Also, let N{v) = J2u lc(u,«)\A“ count the number of bad candidates for which = v. 
Then, using the mass-transport principle. 


E[N{v)] = Y, v) \ A“] = ^ P[C'(u, u) \ = P[5" \ A’'] ^ e. 

u u 

Let O'" denote the event “there exists a good candidate (m, 0) such that Z^ = u”. Using 
the two last estimates we get 


F{D^) ^ P(iV(u) ^ 1) - P(iV(u) ^ 1) ^ 1 - -Sn-e, 


where (5„ —)■ 0 as n —)■ oo. Choosing e small and n large, we have P(D°) > 

On the event there is some non-fixating particle (u, 0) and some inner time t > T 
such that = 0, so vertex 0 is visited by an active particle after time T. Letting 
T —)■ CX), we get P(site 0 not hxating) ^ |, and by the 0-1 law P(0 not hxating) = 1. □ 


|CRS14] For A = cx, Pc ^ 1. 


We start with the following observation about the particle-hole model; 
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The density of holes filled by time t equals the density of particles settled by time t. 

To see that, let A{x,j,y) denote the event that particle {x,j) settles at site y by time 
t, and let m{x,y) = J2j^A{x,j,y)- The density of particles settled by time t is given by 
and the density of holes hlled by time t is given by KJ2y^{y,x). These are 
eqnal by the MTP, proving the identity. 

Proof. Suppose that fi < 1. Using the above identity, 

P[o contains an unhlled hole at time t] ^ 1 — /i > 0. 

Since the above event is decreasing in t and the lower bound is uniform with respect to t, 

P[o is never visited] > 0. 

In particular, the probability that 0 is visited hnitely many times in the particle-hole 
model is positive, so by the 0-1 law, P[mj,(,(0) < cxo] = 1, which means that the ARW 
with A = cx) hxates. □ 


|CRS14] No hxation when p, = 1. 


The proof is done in full details details in |CR,S13j for two-type annihilating systems, 
that is, with reaction A + B ^ tj). Active particles are particles of type A, and holes are 
particles of type B. The reaction happens when a particle settles, thereby hlling a hole. 

Sketch of the proof. Suppose that sites hxate in the particle-hole model. By the hrst 
result presented in this section, site hxation implies that every particle eventually settles. 
Letting t —)■ cx) in the previous observation, 

P[0 is ever visited] = p. 

In the sequel we will show that, under the assumption of site hxation, 

P [0 is never visited] > 0. 

This implies that p < 1, therefore proving that there cannot be hxation for p = 1. 
Assuming site hxation, necessarily, there exists k eN such that 

P[the number of particles which ever visit 0 equals fc] > 0. 
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Moreover, there exist xi,..., G such that ¥\U] > 0, where 

U = [the particles which ever visit 0 are initially at the sites xi,..., x^]. 


Consider two copies ^ and ^ of the system, coupled as follows. We take X = X, and 
fjQ^x) = TjQ^x) for X ^ {xi,... ,Xfc}. Finally, for x G {xi,... ,Xfc}, we sample po and po 
independently. Now notice that 


P 


U occurs for 4, and po(^i) = • ■ ■ = '^o(^fc) = 0 


= P 


U occurs for ^ 


X P 


Vo{xi) = ■■■ = T]o{xk) = 0 for ^ 


> 0 . 


On the above event, no particle ever visits 0 in the system Indeed, on the above 
event, the initial conhguration of ^ is the same as that of ^ except for the deletion of the 
particles present in {xi,..., x^}. In particular, all the particles which visit the origin in ^ 
are deleted in Recalling that ^ and 4 share the same putative trajectories, by following 
how the effect of deleting such particles propagates in the system evolution, we see that 
in the system ^ no particles can ever visit 0. □ 
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